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ABSTRACT 

Aims. The diffuse morphology and transient nature of coronal mass ejections (CMEs) make them difficult to identify and track using 
traditional image processing techniques. We apply multiscale methods to enhance the visibility of the faint CME front. This enables 
an ellipse characterisation to objectively study the changing morphology and kinematics of a sample of events imaged by the Large 
Angle Spectrometric Coronagraph (LASCO) onboard the Solar and Heliospheric Observatory (SOHO) and the Sun Earth Connection 
Coronal and Heliospheric Investigation (SECCHI) onboard the Solar Terrestrial Relations Observatory (STEREO). The accuracy of 
these methods allows us to test the CMEs for non-constant acceleration and expansion. 

Methods. We exploit the multiscale nature of CMEs to extract structure with a multiscale decomposition, akin to a Canny edge 
detector. Spatio-temporal filtering highlights the CME front as it propagates in time. We apply an ellipse parameterisation of the front 
to extract the kinematics (height, velocity, acceleration) and changing morphology (width, orientation). 

Results. The kinematic evolution of the CMEs discussed in this paper have been shown to differ from existing catalogues. These 
catalogues are based upon running-difference techniques that can lead to over-estimating CME heights. Our resulting kinematic 
curves are not well-fitted with the constant acceleration model. It is shown that some events have high acceleration below ~5 Rq. 
Furthermore, we find that the CME angular widths measured by these catalogues are over-estimated, and indeed for some events our 
analysis shows non-constant CME expansion across the plane-of-sky. 

Key words. Sun: coronal mass ejections (CMEs) - Sun: activity - Techniques: image processing - Methods: data analysis. 



1. Introduction 

Coronal mass ejections (CMEs) are large-scale eruptions of 
plasma and magnetic field with energies up to, and be- 
yond, 10 32 ergs. The mass of particles expelled in a CME 
can amount to 10 15 g, typically traveling from the Sun 
at velocities of hundreds up to several thousand kilome- 
ters per second dGoslinget ail 119761; iHundhausen et al.l 1994; 
iGopalswamv et alj|200lb iGallagher et alj 2003b . Thev can lead 
to significant geomagnetic disturbances on Earth and may have 
negative impacts upon space-borne instruments susceptible to 
high levels of radiation when outside the protection of Earth's 
magnetosphere. While they have been investigated by remote 
sensing and in situ measurements for more than 30 years, their 
kinematic and morphological evolution through the corona and 
interplanetary space is still not completely understood. 

It is well known that CMEs are associated with filament 
erupt ions and solar flares dZhang & Wand 120021 ; iMoon et aU 
120041) but the driver mechanis m remains elusive . Several theoret- 
ical models are described by iRrall et aL ( 2001 ) with in the con 



kinematics of an erupting flux-rope can be described using a 
force-balance equation, which includes gas pressure, gravity 



rorce-balance equation, wmcn includ es gas pressure, gr avity 
and the Lorentz force dKrall & Chenl 120051; IChen et al.l 12006; 



Chenl 



1996). These in- 



text of the magnetic flux-rope model of 1 
elude flux injection, magnetic twisting, magnetic energy release 
and hot plasma injection. The three-dimensional magnetic flux- 
rope model, in which a converging flow toward the neutral line 
result s in reconnection b e neath the flux-ro pe dForbes &" Priest 
119951: lAmari et alj 120031; iLin et alj 120071) . assumes that the 
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Manchester 2007). An alternative to this is the magnetic break- 
out model in which the CME eruption is triggered by reconnec- 
tion between the overlying field and a neighbouring flux sys- 
tem dAntiochosetalJll999l iLynch et alJl2004t iMacNeice etall 
2004). The increased rate of outward expansion drives a faster 
rate of breakout reconnection yielding the positive feedback re- 
quired for an explosive eruption. The models are dependent on 
geometrical properties of the CME, such as its width and radius, 
and they are designed to give an indication of the processes that 
drive CME kinematics. Thus it is important to develop methods 
of localising the CME front and characterising it in the observed 
data with high accuracy for model comparisons. 

To date, most CME kinematics are derived from running- 
difference images whereby each image is subtracted from the 
next in order to highlight regions of changing intensity. Once the 
CME is identified, either manually through point-and-click es- 
timates, or by a form of image thresholding and segmentation, 
height-time plots are produced and the velocity and accelera- 
tion of each event determined. CMEs are often thought of as 
being either gradual or impulsive depending on both their erup- 
tion mechanis m (e.g. prominence lift off or flare driven) and 
their speed dSheelev et alj|1999t IMoon et alj|2002l) . The faster 
events tend to have an average negati ve acceleration attributed to 
aerod ynamic drag of the sola r wind jC argill 2004; Vrsna k et al. 
120041) . It was also shown by Gallaghe ret all d2003l) andlZhang 
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d2005l) that CMEs may undergo an early impulsive acceleration 
phase below ~2 Rq. 

Current methods of CME detection have their limitations, 
mostly since these diffuse objects have been difficult to iden- 
tify using traditional image processing techniques. These diffi- 
culties arise from the transient nature of the CMEs, the scat- 
tering effects and non-linear intensity profile of the surround- 
ing corona, and the interference of cosmic rays and solar en- 
ergetic particles (SEPs) appearing as noise on the coronagraph 
detectors. Observations made by the Solar and Heliospheric 
Observatory's (SOHO) Large Angle Spectrometric Coronagraph 
(LASCO) are compiled into a CM E catalog at the Coor dinated 
Data Analysis Workshop (CD AW; lYashiro et aT] |2004|) which 
operates by tracking the CME in LASCO/C2 and C3 running 
difference images to produce height-time plots of each event. It 
is a wholly manual procedure and is subject to user bias in in- 
terpreting the data. The Computer Aided CME Tracking rou- 
tine (CACTus; iRobbrecht & Berghmansl [20041) is based upon 
LASCO/C2 and C3 running difference images. The images are 
unwrapped into polar coordinates and angular slices are stacked 
together in a time-height plot. CMEs thus appear as ridges in 
these plots, detected by a Hough Transform. The nature of this 
detection constrains the CMEs to have constant velocity and 
zero acce leration. The Solar Eruptive Event Detection Systems 
(SEEDS', Olm edo et alj2008b is an automatic detection based on 
LASCO/C2 running difference images unwrapped into polar co- 
ordinates. The algorithm uses a form of threshold segmentation 
to approximate the shape of the CME leading edge, and automat- 
ically determines the height, velocity and acceleration profiles. 

In this work we apply a new multiscale method of analysing 
CMEs. The use of multiscale methods in astr ophysics have 
proven effective at denqising spectra and images (Murtag h et alj 
ll995tlFligge & Solankill 997). analysing solar active region evo- 
lution (jHewett e t al. 2008|), and enhancing solar c oronal images 
dStenborg & Cobellil 120031; IStenborg et all 120081) . Multiscale 
analysis has the benefit of working on independent images 
without any need for differencing. A parti cular application of 
multis cale decompositions, implemented by Young & Gallagher 
(2008), uses high and low pass filters convolved with the im- 
age data to exploit the multiscale nature of the CME. This high- 
lights its intensity against the background corona as it propa- 
gates through the field-of-view, whil e neglecting small scale fea- 
tures (essentially denoising the data). Young & Gallagher (2008) 
also describe the use of nonmaxima suppression to trace the 
edges in the CME images, and they show the power of multi- 
scale methods over previous edge detectors such as Roberts and 
Sobel. With these methods for defining the front of the CME we 
can characterise its morphology (width, orientation) and kine- 
matics (position, velocity, acceleration) in coronagraph images. 
Follo wing previous discussions on the errors in CME heights 
(e.g. IWenet ail 120071) . multiscale methods allows us to deter- 
mine the kinematics to a high degree of accuracy in order to 
confidently compare our results to theory. While these methods 
are not currently automated, they have great potential for such 
an implementation which the authors intend to pursue. 

In Sect. 2 we discuss current image pre -processing me thods 
from the c oronagraphs onboa rd the SOHO (Domingo 1997b and 
STEREO dKaiser et al.ll2008l) spacecraft. We outline the imple- 
mentation of multiscale methods and ellipse fitting to charac- 
terise the CME front. In Sect. 3 we present a sample of CME 
events analysed by these methods, and Sect. 4 is a discussion of 
these first results and conclusions. 




Fig. 1. Raw (left) and pre-processed image (right) of a CME ob- 
served by LASCO on 2004 April 1 . The pre-processing includes 
median filtering, background subtraction and occulter masking. 



2. Observations & Data Analysis 

In this paper a sample of CMEs observed by SOHO/LASCO 
and STEREO/SECCHI are studied, namely the gradual events 
of 2000 January 2, 2000 April 18, 2001 April 23, 2004 April 
1, 2007 October 8 and 2007 November 16, and the impulsive 
events of 2000 April 23 and 2002 April 21. 

The LASCO/C2 and /C3 coronagraphs dBrueckner et alj 

fl995l) onboard SOHO give a nested field-of-view extending 
from approximately 2. 2 - 30 Rp). The CO R1/2 coronagraphs 
of the SECCHI suite dHoward et all 12000) onboard STEREO 
image the corona from 1.4 - 15 R0. The image quality of the 
coronagraph observations can be diminished for many reasons, 
including instrumental effects (e.g. scattered light), noise from 
cosmic rays and SEPs, or data dropouts. In order to minimise 
these effects, we use the following standard preprocessing meth- 
ods. Firstly the images are normalized with regard to exposure 
times in order to correct for temporal variations in the image 
statistics. Secondly a median filter is applied to remove pixel 
noise, replacing hot pixels with a median value of the surround- 
ing pixel intensities. Finally we perform a background subtrac- 
tion, obtained from the minimum of the daily median images 
across a time span of a month, and remove the occulting disc 
with a zero mask. These steps lead to a clear improvement in the 
image quality for CME study (Fig.[T|), after which we apply our 
methods of multiscale analysis. 

In recent years the use of wavelets has be en increasingly evi- 
dent in image processing of solar structures (|Stenborg & Cobellil 
120031; IStenborg et alj|2008t lYoung & Gallagherj|2008l) . Here we 
illustrate briefly the formalism of the ID wavelet transform, 
and discuss how this is extended to 2D. (For a more detailed 
discussion on wav elet transforms the re ader is directed to e.g. 
iBurrus et all d!998l) and lLin & Liul d2004l) .) The remainder of the 
section outlines our methods of using the multiscale analysis to 
define the CME front for characterisation. 

2. 1 . Multiscale Filtering 

The continuous wavelet transform (CWT) of a signal f(x) e 
L 2 (R) with respect to the mother wavelet if/(x) e L 2 (R) is de- 
fined by: 

W(a,b) = J f(x)if, a , b (x)dx (1) 

where ^/ a ,b(x) is a spatially localised function given by: 
1 x-b 

*ajb&) = -plK ) (2) 
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Fig. 2. Top left, the horizontal, and top right, the vertical coef- 
ficients from the high-pass filtering at scale 3. Bottom left, the 
corresponding magnitude (edge strength) and bottom right, the 
angle information (0 - 360°) taken from the gradient space, for 
a CME observed in LASCO/C2 on 2004 April 1. 
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(b) 01:00 UT 



Fig. 3. The vectors plotted represent the magnitude and angle 
determined from the gradient space of the high-pass filtering at 
scale 3. The CME of 2004 April 1 shown here is highlighted 
very effectively by this method. 



and a, b e L 2 (R) correspond to the scaling (dilation) and shifting 
(translation) of i/^x). Note that the CWT spectrum obtained has 
no restriction as to how many scales are used, nor of the spacing 
between the scales. 

On the other hand, the discrete wavelet transform (DWT) 
consists of a scaling function <p(x) and corresponding wavelet 
if/(x) with finite support [0, /] {I being a positive number) given 
by: 



4>{x) = V2 Y^h s cf>(2x- s) 



(3) 



iff{x) = V2V^(2x-s) 



(4) 



s=0 



where h and g are constants called the low-pass and high-pass 
filter coefficients respectively. Unlike the CWT, the DWT de- 
composes a signal into a set of wavelet components with dyadic 
scaling, i.e. the scale increases in powers of 2 (2 1 , 2 2 , etc). In this 
way, the computational time required to perform the decompo- 
sition is drastically reduced. 

We explore a method of multiscale decomposition by means 
of the DWT in 2D through the use of low and high pass filters 
using a discrete approximation of a Gaussian 6 and its derivative 
ift respectively. Since 8(x, y) is separable we can write: 



dx 



^ y (x,y) = 9(x) 



d6(y) 
dy • 



(5) 



(6) 



Successive convolutions of an image with the filters produces the 
scales of decomposition, with the high-pass filtering providing 
the wavelet transform of image I(x, y) in each direction: 



WJ = WJ{x,y) = iff x (x,y)*I(x,y) 
W y I = W y I(x,y) = t/r y (x,y)*I{x,y). 



(7) 
(8) 



Akin to a Canny edge detector ( Young & Gallagher 2008), these 
horizontal and vertical wavelet coefficients are combined to form 
the gradient space F for each scale: 



T(x,y) = [WJ, W y l]. 



(9) 



The gradient information has an angular component a and a 
magnitude (edge strength) M: 



a(x,y) = tan' 1 (w,I / Wj) 



M(x, y) 



+ (Wyl) 2 



(10) 



(ID 



s=0 



The resultant horizontal and vertical detail coefficients, and the 
magnitude and angular information are illustrated in Fig. [2] 

Overlaying a mesh of vector arrows on the data shows how 
the magnitude and angular information illustrate the progression 
of the CME. Each vector is rooted on a pixel in the gradient 
space, and has a length corresponding to the magnitude M with 
an angle from the normal a (Pig. [3j. By utilising all the informa- 
tion gained from the convolutions on the data array it becomes 
possible to threshold out the CME with a view to characterising 
its progression through space. 

The magnitude information was found to have the highest 
signal-to-noise ratio at the third scale of the decomposition. This 
scale is very effective at smoothing unwanted artefacts such as 
cosmic rays which the median filter may have missed. The an- 
gular component a of the gradient specifies a direction which 
points across the greatest intensity change in the data (an edge). 
A threshold is specified with regard to this gradient direction in 
order to chain pixels along maxima, highlighting the edges in the 
image. A spatio-temporal threshold is implemented with regard 
to changes of magnitude in order to highlight moving features. 
It works by creating a specific detection mask which is used to 
pull out the edges along the CME front. In cases of faint CMEs 
or strong streamer deflections the filter is presently limited by 
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working on only one scale, meaning current analysis involves 
a user removing/including certain edges that the algorithm has 
mistakenly retained/discarded. Future work on more than one 
scale may help alleviate this issue, since the additive angular 
distribution along a CME front is significantly larger than along 
a streamer when considered through multiple scales. Including 
this caveat will better enable the algorithm to automatically de- 
termine the CME edges, leading to a near real-time characteri- 
sation and kinematic analysis. 

2.2. Characterising the CME Front 

Using a model such as an ellipse to characterise the CME front 
across a sequence of images, has the benefit of providing the 
kinematics and morphology of a moving and/or expanding struc- 
ture. The ellipse's multiple parameters, namely its changeable 
axes lengths and tilts, is adequat e for approximati ng the vary- 
ing curved structures of CMEs. IChen et all d 19971) suggest an 
ellipse to be the two-d imensional projection of a flux rope, and 
iKrall & St. Cvrl d2006l) use ellipses to parameterise CMEs and 
explore their geometrical properties. We fit ellipses to the points 
determined to be along the CME front by considering a radial fan 
from Sun centre across the defined edges. This means there are 
more points along the front than on the flanks of the CME for in- 
clusion in the fit, and the edges never double back on themselves 
in cases where the CME's internal structure might be otherwise 
included. A kinematic analysis provides height, velocity and ac- 
celeration profiles; while the ellipse's changing morphology pro- 
vides the inclination angle and angular width. Measuring these 
properties in the observed data is vitally important for accurate 
comparison with theoretical models. 

The implementation of the ellipse fitting routine is based 
upon an initial guess of ellipse centre as the average of the points 
specified along the front. The ellipse equation (in polar coordi- 
nates) is defined as: 



p 2 cos 2 (i) p 2 sin 2 a) 

h = 1 

a 2 b 2 



(12) 



where a and b are the lengths of the semimajor and semiminor 
axes respectively, so allowing for an inclination angle y on the 
ellipse gives: 



2u2 



a l b 



cos(2 W '-2 r ) 



(13) 



where a>' = to + y, as illustrated in Fig. [4] This gives a first ap- 
proximation which can then be used to iteratively float the ellipse 
parameters until a Levenberg-Marquardt least-squares minimi- 
sation is reached. 

CME height measurements are taken as the height of the fur- 
thest point on the ellipse from Sun centre. The angular width is 
taken as the cone angle of the ellipse from Sun centre, and the 
tilt of the ellipse is given by the calculated angle y. Note that 
in cases where the code produces an extremely large and oblate 
ellipse with one apex approximating the CME front, the width 
and tilt information is deemed redundant. Hence the resulting 
analysis of some events can have less data points included in the 
width and tilt plots than in the height-time plots. 

2.3. Error Analysis & Model Fitting 

The front of the CME is determined through the multiscale de- 
composition and consequent rendering of a gradient magnitude 




Fig. 4. Ellipse inclined at angle y, with semimajor axis a, 
semiminor axis b, and radial line p inclined at angle to to the 
semimajor axis. 



space. At scale 3 of the decomposition the smoothing filter is 2 3 
pixels wide, which we use as our error estimate in edge position. 
This error is input to the ellipse fitting algorithm for weighting 
the ellipse parameters, and a final error output is produced for 
each ellipse fit. In the case of a fading leading edge the reduced 
amount of points along the front will increase the error on our 
analysis. The final errors are displayed in the height-time plots 
of the CMEs, and are used in the velocity and acceleration cal- 
culations. The derivative is a 3-point Lagrangian interpolation, 
so there is an enhancement of error at the edges of the data sets. 

The errors on the heights are used to constrain the best fit to 
a constant acceleration model of the form: 



h(t) = aot + vot + ho 



(14) 



where t is time and oq, vo and are the acceleration, initial ve- 
locity and initial height respectively. This provides a linear fit to 
the derived velocity points and a constant fit to the acceleration. 

An important point to note is the small time error (taken to be 
the image exposure time of the coronagraph data) since the anal- 
ysis is performed upon the observed data frames individually. 
Previous methods of temporal-differencing would increase this 
time error. With these more accurate measurements we are bet- 
ter able to determine the velocity and acceleration errors, leading 
to improved constraints upon the data and providing greater con- 
fidence in comparing to theoretical models. 

3. Results 

This section outlines events which have been analyzed using our 
multiscale methods. We use data from the LASCO/C2 and C3, 
and SECCHI/COR1 and COR2 instruments, and preprocess the 
images as discussed in Sect. 2. The ellipse fitting algorithm ap- 
plied to each event gives consistent heights of the CME front 
measured from Sun centre to the maximum height on the el- 
lipse, and these lead to velocity and acceleration profiles of our 
events. The ellipse fitting also provides the angular widths and 
orientations, as shown below. The velocity, acceleration and an- 
gular width results of each method are highlighted in Tables 1, 
2 and 3. In each instance we include the values from CACTus, 
CDAW and SEEDS, noting that CACTus lists a median speed 
of the CME; CDAW provide the speed at the final height and 
from the velocity profile we infer the speed at the initial height; 
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and the SEEDS detection applies only to the LASCO/C2 field- 
of-view but doesn't currently provide a velocity range or profile. 
Note also that the CMEs of 2007 October 8 and 2007 November 
16 are analysed in SECCHI images by CACTus and our multi- 
scale methods (marked by asterisks in the Tables), while CDAW 
and SEEDS currently only provide LASCO analysis. It is clear 
that many of the CACTus, CDAW and SEEDS results lie outside 
the results and error ranges of our analysis. 



Table 1. Summary of CME velocities as measured by CACTus, 
CDAW, SEEDS and multiscale methods. 



Date 


CACTus 


CDAW 


SEEDS 


Multiscale 






km s" 1 






2000 Jan 02 


512 


370 - 794 


396 


396 - 725 


2000 Apr 18 


463 


410-923 


339 


324 - 1049 


2000 Apr 23 


1041 


1490 - 898 


595 


1131 - 1083 


2001 Apr 23 


459 


540-519 


501 


581 -466 


2002 Apr 21 


1103 


2400 - 2388 


702 


2195 -2412 


2004 Apr 01 


487 


300-613 


319 


415-570 


2007 Oct 08 


235* 


85-331 


103 


71 -330* 


2007 Nov 16 


337* 


210-437 


154 


131 -483* 



*analysis of SECCHI data rather than LASCO. 



Table 2. Summary of CME accelerations as measured by 
CACTus, CDAW, SEEDS and multiscale methods. 



Date 


CACTus 


CDAW 

-2 

m s 


SEEDS 


Multiscale 


2000 Jan 02 





21.3 


-5.8 


14.7 


±3.6 


2000 Apr 18 





23.1 


17.5 


32.3 


± 3.5 


2000 Apr 23 





-48.5 


-8.9 


-4.8 


±20.6 


2001 Apr 23 





-0.7 


-1.4 


-4.8 


±4.1 


2002 Apr 21 





-1.4 


33.5 


32.5 


±26.6 


2004 Apr 01 





7.1 


12.9 


4.4 


±2.0 


2007 Oct 08 


0* 


3.4 


2.4 


5.7 4 


:0.9* 


2007 Nov 16 





4.9 


11.0 


13.7 


± 1.7* 



*analysis of SECCHI data rather than LASCO. 



Table 3. Summary of CME angular widths as measured by 
CACTus, CDAW, SEEDS and multiscale methods. 



Date 


CACTus 


CDAW 
degrees 


SEEDS 


Multiscale 


2000 Jan 02 


160 


107 


96 


50 


-95 


2000 Apr 18 


106 


105 


108 


68- 


110 


2000 Apr 23 


352 


360 


130 


96- 


130 


2001 Apr 23 


124 


91 


74 


55 


-60 


2002 Apr 21 


352 


360 


186 


53 


-65 


2004 Apr 01 


66 


79 


58 


44 


-38 


2007 Oct 08 


52* 


82 


59 


23- 


60* 


2007 Nov 16 


68* 


78 


54 


40- 


55* 



*analysis of SECCHI data rather than LASCO. 
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Start Time (02-Jan-00 05:54:18) 

Fig. 6. Kinematic and morphological profiles for the ellipse fit to 
the multiscale edge detection of the 2000 January 2 CME ob- 
served by LASCO/C2 and C3. The plots from top to bottom 
are height, velocity angular width, and ellipse tilt. The CDAW 
heights are over-plotted with a dashed line. The height and ve- 
locity fits are based upon the constant acceleration model. 



3.1. Arcade Eruption: 2000 January 2 

This CME was first observed in the south-west at 06:06 UT on 
2000 January 2 and appears to be a far-side event associated with 
an arcade eruption consisting of one or more bright loops. 

The height-time plot has a trend not unlike that of CDAW 
(overplotted in top Fig.|6]with a dashed line). However the offset 
of the CDAW heights - which puts them outside our error bounds 
- may be due to how the difference images are scaled for display. 
This is a problem multiscale methods avoid. From Fig. |6l the 
velocity-fit was found to be increasing from 396 to 725 km s , 
giving an acceleration of 14.7 ± 3.6 m s~ 2 . The ellipse fit spans 
approximately 50 - 70° of the field-of-view in the inner portion 
of C2, and expands to over 95° in C3. This expansion may sim- 
ply be attributed to the inclusion of one or more loops in the 
ellipse fit as the arcade traverses the LASCO/C2 and C3 fields- 
of-view. The orientation of the ellipse as a function of time is 
shown in bottom Fig. [6] It can be seen that the orientation angle 
of the CME increases to approximately 100° before decreasing 
toward 60°. 

The constant acceleration model is not a sufficient fit to the 
data in this event. The kinematics produced from the multiscale 
edge detection would be better fit with a non-linear velocity and 
a non-constant acceleration. This would show the CME to have 
a period of decreasing acceleration in the C2 field-of-view, lev- 
eling off to zero in C3 (if not decelerating further). 
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2000 January 2 2000 April 18 2000 April 23 2001 April 23 




Fig. 5. A sample of ellipse fits to the multiscale edge detection of the events studied. For each event the upper and lower image show 
LASCO/C2 and /C3 (except for the 2007 events which show SECCHI/COR1 and /COR2). 
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18:00 19:00 30:00 21:00 
Start Time (18-Apr-00 15:47:37) 



Fig. 7. Same as Fig.|6]but for the gradual CME of 2000 April 18. 



3.2. Gradual/Expanding CME: 2000 April 18 

This CME was first observed off the south limb at 16:06 UT on 
2000 April 18 and exhibits a flux rope type structure. 

The height-time plot for this event has a trend similar to that 
of CDAW (overplotted in top Fig. [7] with a dashed line). The 
velocity-fit was found to be linearly increasing from 324 to over 
1 000 km s _1 , giving an acceleration of 32.3 ± 3.5 m s~ 2 . The 
ellipse fit spans from 68° of the field-of-view in the inner portion 
of C2, to approximately 1 10° in C3. The orientation of the ellipse 
as a function of time is shown to increase from just above 0° to 
over 60° in Fig. [7J 

This CME is fitted well with the constant acceleration model 
but shows an increasing angular width implying expansion 
across the field-of-view. 



3.3. Impulsive CME: 2000 April 23 

This impulsive CME was first observed in the west at 12:54 UT 
on 2000 April 23 and exhibits strong streamer deflection. 

The height-time plot derived using our methods has a trend 
which diverges from that of CDAW (overplotted in top Fig. [8] 
with a dashed line). The velocity-fit was found to be linearly de- 
creasing from 1 131 to 1 083 km s _1 , giving a constant decelera- 
tion of -4.8 ± 20.6 m s~ 2 . The CME is present for one frame in 
C2 with an ellipse fit spanning 96°, increasing to approximately 
120 - 130° in the C3 field-of-view, and the orientation of the el- 
lipse as a function of time is shown to rise from 71° to 95° then 
fall to 64° (see bottom Fig. [8]). 

This event is modeled satisfactorily with a constant deceler- 
ation. However, due to the impulsive nature of the CME there 
are only a few frames available for analysis, making it difficult 
to constrain the kinematics. 
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Fig. 8. Same as Fig. |6]but for the impulsive CME of 2000 April 
23. 



3.4. Faint CME: 2001 April 23 

This CME was first observed in the south-west at 12:39 UT on 
2001 April 23 and exhibits some degree of streamer deflection. 

The height-time plot has a similar trend to CDAW (overplot- 
ted in top Fig. [9] with a dashed line). The velocity-fit was found 
to be linearly decreasing from 58 1 to 466 km s , giving a decel- 
eration of -4.8 ±4.1 m s~ 2 . The ellipse fit spans approximately 
55 - 60° of the field-of-view throughout the event, and the ori- 
entation of the ellipse as a function of time is shown to decrease 
from approximately 50° to almost 0° (see Fig.[9]l. 

This CME is fitted well with the constant acceleration model. 



3.5. Fast CME: 2002 April 21 

This CME was first observed in the west from 01 :27 UT on 2002 
April 21. 

The height-time plot follows a similar trend to that of CDAW 
(overplotted in top Fig. [10] with a dashed line). The velocity -fit 
was found to be linearly increasing from 2 195 to 2412 km s , 
giving a constant acceleration fit of 32.5 ± 26.6 m s~ 2 . The el- 
lipse fit spans 53° in C2, and shows an likely increasing trend 
to 65° in C3. The orientation of the ellipse as a function of time 
is shown to scatter about 1 15° though it drops to approximately 
8 1 in the final C3 image. 

The kinematics of this event are not modeled satisfactorily 
by the constant acceleration model, since the fits do not lie within 
all error bars. The argument for a non-linear velocity profile, 
with a possible early decreasing acceleration, is justified for this 
event, although the instrument cadence limits th e data set avail- 
able for interpretation. The previous analysis of Gallagher et al. 
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Fig. 9. Same as Fig. ©but for the faint CME of 2001 April 23. Fig. 10. Same as Fig. ©but for the fast CME of 2002 April 21. 



d2003l) resulted in a velocity of -2 500 km s" 1 past -3.4 R© 
which is consistent with our results past ~6 R© in Fig.fTUl 

3.6. Flux-Rope/Slow CME: 2004 April 1 

This CME was first observed in the north-east from approxi- 
mately 23:00 UT on 2004 April 1, is in the field-of-view for 
over 9 hours, and exhibits a bright loop front, cavity and twisted 
core. 

The height-time plot follows a similar trend to that of CDAW 
(overplotted in top Fig. QT|with a dashed line). The velocity-fit 
was found to be linearly increasing from 415 to 570 km s , giv- 
ing an acceleration of 4.4 ± 2.0 m s~ 2 . Note also that the kine- 
matics of this event exhibit non-linear structure clearly seen in 
the velocity and acceleration profiles. The ellipse fit spans ap- 
proximately 44° in C2, stepping down to approximately 38° in 
C3. The orientation of the ellipse as a function of time is shown 
to jump down from approximately 130° in C2 to approximately 
70-80° in C3. 

This event shows unexpected structure in the velocity and 
acceleration profiles which indicates a complex eruption not sat- 
isfactorily modeled with constant acceleration. 

3. 7. STEREO-B Event: 2007 October 8 

This CME was first observed in the west from approximately 
12:00 UT on 2007 October 8, and is best viewed from the 
STEREO-B spacecraft. It is noted that the kinematics as mea- 
sured by S OHO and ST EREO will be different due to projec- 
tion effects ( Vrsnak et al. 20Q7|). It is intended that these effects 
be measured and corrected for in future work, especially as the 



spacecraft separations increase over time. On this date STEREO- 
B was at an angular separation of 16.5° from Earth. 

The height-time plot for this event is shown in Fig. Q~2] 
The velocity-fit was found to be linearly increasing from 71 to 
330 km s , giving an acceleration of 5.7 + 0.9 m s~ 2 . The el- 
lipse fit in COR1 spans approximately 23° stepping up to a scat- 
ter about 40 - 50° which rises slightly to 50 - 60° in COR2. 
The orientation of the ellipse as a function of time is shown to 
increase from 55 - 110° then jumps to an approximately steady 
scatter about 180- 190°. 

This CME is fitted well with the constant acceleration model. 

3.8. STEREO-A Event: 2007 November 16 

This CME was first observed in the west from approximately 
08:26 UT on 2007 November 16, and is best viewed from the 
STEREO-A spacecraft. On this date STEREO-A was at an an- 
gular separation of 20.3° from Earth. 

The height- time plot for this event is shown in Fig. [13] The 
velocity-fit was found to be linearly increasing from 131 to 
483 km s _1 , giving an acceleration of 13.7 ± 1.7 m s~ 2 . The 
ellipse fit in COR1 spans approximately 40 - 50° stepping up 
slightly to a scatter about 45 - 55° in COR2. The orientation of 
the ellipse as a function of time is shown to start at 153° and end 
at 120° with the mid points scattered about 170°. 

This CME is fitted well with the constant acceleration model. 



4. Discussion & Conclusions 

A variety of theoretical models have been proposed to describe 
CMEs, especially their early propagation phase. Observational 
studies, such as those outlined above, are necessary to determine 
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Fig. 11. Same as Fig.|6]but for the slow CME of 2004 April 1. Fig. 13. Same as Fig.[[2]but for the CME of 2007 November 16. 
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Fig. 12. Same as Fig. but for the CME of 2007 October 8 ob- 
served by SECCHI/COR1 and COR2. For this reason CDAW is 
not overplotted. 



CME characteristics. We argue that the results of previous meth- 
ods are limited in this regard due mainly to large kinematic errors 
which fail to constrain a model, an artefact of CME detection 
based upon either running- (or fixed-) difference techniques or 
other operations. Current methods fit either a linear model to the 
height-time curve, implying constant velocity and zero acceler- 
ation (e.g. CACTus) or a second order polynomial, producing a 
linear velocity and constant acceleration (e.g. CDAW, SEEDS). 
The implementation of a multiscale decomposition provides a 
time error on the scale of seconds (the exposure time of the in- 
strument) and a resulting height error on the order of a few pix- 
els. The height- time error is used to determine the errors of the 
velocity and acceleration profiles of the CMEs. It was shown that 
for certain events the results of CACTus, CDAW and SEEDS can 
differ significantly from our methods, as illustrated in the Tables 
of Sect. 3. 

Our results clearly confirm that the constant acceleration 
model may not always be appropriate. The 2000 January 2 and 
2002 April 21 CMEs are good examples of the possible non- 
linear velocity profile and consequent non-constant acceleration 
profile (see Fig. [6] and Fig. [TOt. Indeed these events are shown 
to have a decreasing acceleration, possibly to zero or below, as 
the CMEs traverse the field-of-view. Simu lations of the break- 
out model outlined in iLvnch et"aT] d2004l) resulted in constant 
acceleration fits which do not agree with these observations. It 
may be further noted that the events of 2001 April 23 and 2004 
April 1 show a possible decreasing acceleration phase early on, 
though within errors this cannot be certain (see bottom Fig.|9]and 
Fig. ITTb . Furthermore, the structure seen in some events would 
indicate that the CME does not progress smoothly. The veloci- 
ties of the 2004 April 1 CME in Fig. HUand the 2007 November 
16 CME in Fig. [13] show non-smooth profiles and may imply 
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a form of bursty reconnection or other staggered energy release 
driving the CME. Other profiles such as Fig. [6] and to a lesser 
extent Figs. [9] and [10] may show a stepwise pattern, indicative 
of separate regimes of CME progression. None of the current 
CME models indicate a form of non-smooth progression, al- 
though the flux-rope model does describe an early acceleration 
re gime giving a no n-linear velocity to the eruption (see Fig. 1 1.5 
in lPriest & Forbesll2000l) . 

It may be concluded that the angular widths of the events are 
indicative of whether the CME expands radially or otherwise 
in the plane-of-sky. For the CMEs studied above, the observa- 
tions of 2000 April 18, 2000 April 23 and 2002 April 21 show 
a super-radial expansion (see Fig. [71 Fig. [8] and Fig. ITOb. These 
events also show high velocities, obtaining top speeds of up to 
1 000 km s _1 , over 1 100 km s _I and 2500 km s _1 respectively, 
and may therefore indicate a link between th e CME expan- 
sion a nd speed. Furthermore, it is suggested bv lKrall&St.Cvrl 
(2006) that the flux-rope model can account for different ob- 
served expansion rates due to the axial versus broadside view 
of the erupting flux system. 

The observed morphology of the ellipse fits may be further 
interpreted through the tilt angles plotted in Sect. 3. In knowing 
the ellipse tilt and the direction of propagation of the CME it 
is possible to describe the curvature of the front. For the events 
above, the changing tilt and hence curvature is possibly signif- 
icant for the 2000 April 18, 2004 April 1 and 2007 October 8 
events (see bottom Fig.|7J Fig.[TT]and Fig.[T2ii. The elliptical flux 
rope model of Kra ll et al.l d2006h was shown to have a chang- 
ing orientation of the magnetic axis which results in a dynamic 
radius of curvature of the CME, possibly accounting for these 
observed ellipse tilts. 

The work outlined here is an initial indication that the zero 
and constant acceleration models in CME analysis are not an 
accurate representation of all events, and the over-estimated an- 
gular widths are not indicative of the true CME expansion. The 
ellipse characterisation has provided additional information on 
the system through its changing width and orientation which 
shall be further investigated with the analysis of many more 
events. This work will be further explored and developed with 
STEREO data whereby the combined view -points can give addi- 
tional kinematic constraints and may lead to a correction for pro- 
jection effects and possible three dimensional reconstructions. 
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